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Abstract 

We investigate the disorder-induced localization transition in Bose-Einstein condensates for the 
Anderson and Aubry-Andre models in the non-interacting limit using exact diagonalization. We 
show that, in addition to the standard superfluid fraction, other tools such as the entanglement 
and fidelity can provide clear signatures of the transition. Interestingly, the fidelity exhibits good 
sensitivity even for small lattices. Effects of the system size on these quantities are analyzed in 
detail, including the determination of a finite-size-scaling law for the critical disorder strength in 
the case of the Anderson model. 

PACS numbers: 03.75.Nt, 67.85.Hj, 64.60.an, 72.15.Rn 
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I. INTRODUCTION 



Localization of waves in a disordered medium was first predicted by Anderson [f ] in the 
context of non-interacting electrons in a random crystal. Recently, developments in the 
research into ultracold atomic gases have enlarged the possibilities for studying Anderson 
localization in new systems. A breakthrough in the area was the experimental realization 
of Anderson localization of a Bose-Einstein condensate (BEC) in two different kinds of 
disordered optiea. potentials QQ: Bi.ly et a,. 3 employed a BEC of -Rb atoms in the 
presence of a controlled disorder produced by a laser speckle, while Roati et al. [3J utilized a 
BEC of essentially non-interacting 39 K atoms in combination with a one- dimensional quasi- 
periodic lattice to observe Anderson localization. 

Subsequently, this issue has been attracting increasing interest from both the experi- 
mental and theoretical communities. From the experimental point of view, a variety of 
techniques have been implemented, such as bichromatic optical lattices 13-51, speckle laser 

flflQ ft 

patterns |2j, la, 17J , and disordered cold atom lattices |8| . Theoretically, different approaches 



such clS db variational method 
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9], quantum Monte Carlo simulation [10|, 111], exact diago- 
Q,Bfl, renormalization^group QQ, density-matrix renormalization group 

23], and a perturbative treatment (24 1, among others, 



ization 

mean-field approximations 
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have been explored in this context of disordered ultracold atoms. 

The usual model for bosons on a lattice utilizes the so-called Bose-Hubbard Hamiltonian, 
with the general form 



H = ti^i + VtY^{a\a-j + OjCLi) + ^ ^ n^rii - 1) 

i (ij) i 



(1) 



with the standard notation for creation, annihilation, and number operators for bosons at 
lattice sites. Each site is considered to have a single bound state of energy the hopping 
between sites is restricted to nearest neighbors, with amplitude Q, and U is a local repulsive 
interaction. 

Randomness is introduced via the values of £j. These values can have a truly random 
(usually uniform) distribution in the range —A/2 < Si < A/2, which characterizes the 
Anderson model jl|. Alternatively, they can vary periodically with a period incommensurate 
with the lattice spacing. This is the case in the Aubry- Andre (AA) model 25] for which the 
energies, in the one- dimensional written as Si = A cos(27r/3i), where f3 = (l + y / 5)/2 
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is the golden ratio, and % assumes integer values from 1 to L (the system size in units of the 
lattice spacing). Experimentally, speckle patterns can be modeled by the Anderson model, 
while quasi-periodic potentials can be described by the AA model. 

The phenomenology of model ([T]) at low temperatures involves Bose-Einstein condensa- 
tion, superfluidity, a Mott (sapped) phase, and Anderson localization, depending on the 



relative values of f2, U, and A 



261 ] . Theoretical studies of such a rich phenomenology are 



faced with the usual challenges of interacting many-body problems. For arbitrary interac- 
tion strengths perturbation theory cannot be used, which more or less restricts the approach 
to numerical solutions on finite-size lattices. Then, finite-size effects tend to smear out the 
sharp phase boundaries expected in the thermodynamic limit. This is true even for the 
pure Anderson transition in the non-interacting limit. For instance, in his original paper l| 
Anderson proved that his model presents localization for any A > in one dimension, while 
Aubry and Andre 25] proved that in their model localization only occurs for A > 2, but 
these critical values are not clearly seen in small lattices. 

Our aim in this paper is to discuss in detail the effects of finite size in the non-interacting 
limit for a one-dimensional lattice. We focus on the size dependence of quantities that pro- 
vide signatures of quantum phase transitions (QPTs), here applied to Anderson localization. 
In addition to the superfluid fraction, common in studies of Bose-Einstein condensates, we 
also employ quantities borrowed from quantum information theory, like entanglement and 
fidelity 28[]. All of these quantities allow us to obtain the correct thermodynamic limit of 
the critical disorder strength for both Anderson and Aubry-Andre models. In particular, we 
are able to show that the critical disorder strength in the Anderson model obeys a scaling 
law with the system size. 



II. LOCALIZATION TRANSITION 



We consider one-dimensional lattices of L sites, and perform exact numerical diagonal- 
ization of the Hamiltonian matrix, obtaining a chosen number of lowest-energy eigenvalues, 
and the corresponding eigenvectors. We will focus only on ground-state results. Since in- 
teraction is not taken into account here, we restrict our numerical calculations to a single 
particle. When evaluating the quantities of interest for the Anderson model, it is necessary 
to average over a large number of random configurations of the local energies in order to 
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have physically meaningful results. 

For each lattice size, we calculate the superfluid fraction, ground-state entanglement and 
fidelity, as discussed in detail below. This is done for a wide range of disorder strength 
values, with the aim of determining a critical value above which the states are localized. 
Throughout the paper, energy quantities (like A) are measured in units of the tunneling 
amplitude Q, while the lattice size L is expressed in units of the lattice spacing, which is 
equivalent to the number of sites. 



A. Superfluid fraction 

The ground state of a uniform Bose-Einstein condensate should be superfluid. There is 
some controversy as to whether this remains true in the non- interacting limit. As discussed 



by Lieb et al. [29| , if one defines a superfluid as a fluid that does not respond to an external 
velocity field (e.g., rotation of the container walls), then the ideal Bose gas is a perfect 
superfluid in its ground state. This is the approach we adopt here. On the other hand, 
the presence of disorder tends to alter the nature of the quantum states, which become no 
longer extended over the entire system, but localized over a finite distance. Then, the so 
called superfluid fraction (here denoted by / s ) can signal whether the states are extended 
(f s ^0) or not (/, = 0). 

The superfluid fraction for a discrete lattice system is proportional to the difference in 
ground-state energy between the systems with periodic and with twisted boundary condi- 



tions [4] . Instead, we may keep periodic boundary conditions and use a twisted Hamiltonian 

L 

H e = Yl £ i n i + n XX a l a i e ~^ + aj a * e ^) > ( 2 ) 

where 9 is the twist angle. Ground-state energies of the Hamiltonian (J2J) are calculated for 
9 = (E ) as well as for finite 9 (Eg), and the superfluid fraction is given by jj] 

L 2 (Eg - go) 

Js Ntt 9 2 ' 1 } 

being independent of the twist angle 9 as long as 9 « ir. In practice, f s remains the same 
even for twist angles that are a significant fraction of n. Typically, we use 9 = 0.1 in our 
numerical calculations. 
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FIG. 1: (Color online) Average superfluid fraction in the Anderson model as a function of the 
disorder strength A for different lattice sizes. The averages were calculated from 5000 random 
configurations of the potential. The lines correspond to fittings with Eq. @. 

In Fig. [T]we show the superfluid fraction for the Anderson model as a function of disorder 
strength for different lattice sizes. When A = we always have f s = 1, i.e., the system is in 
a superfluid phase. For large enough A we have f s — > 0, and the ground-state is localized. 

We define a characteristic value A^ of the disorder strength for each lattice size by 
observing that the superfluid-fraction curves are well fitted by the function 

f s ~ e ( "^ r . (4) 

For all lattice sizes L we find that the fitting parameter a has approximately the same value 
a ~ g. Lines in Figure [1] show the fittings with Eq. fll]) while symbols are the corresponding 
values obtained by the numerical solution. Fig. [1] seems to indicate that A^ — > when 
L — > oo, as expected. This behavior will be verified in Sec. II III through a finite-size-scaling 
analysis. 

Figure [2] shows the superfluid fraction for the Aubry- Andre model as a function of the 
"disorder strength" A for different lattice sizes. Actually, A is the amplitude of the Aubry- 
Andre potential, i.e., the incommensurate periodic modulation of local energies. We can see 
that f s {A) = for A larger than a characteristic Al, with Al m 2, consistent with the 
expected critical value A c = 2. However, this happens for certain lattice sizes, but for other 
values of L the behavior of f s with A is completely different, even for very similar sizes, as 
exemplified for L = 54 and 55 in Fig. [2j The reason for this apparently puzzling feature 
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FIG. 2: (Color online) Superfluid fraction in the AA model as a function of the potential amplitude 
A for chosen lattice sizes. The expected critical value A c = 2 is well verified for suitably chosen 
lattice sizes, but not for others (e.g., L = 54 in this figure), as discussed in the text. The lines 
indicated at the top right corner are essentially indistinguishable, except near A = 2, where 
rounding is more noticeable for smaller sizes. 

lies in the mismatching between periodic boundary conditions with integer period L and the 
Aubry- Andre potential with irrational period One can verify that Al ~ 2 occurs when 
L belongs to the sequence of Fibonacci numbers |27j. The first members of this sequence, 
for L > 8, are 8, 13, 21, 34, 55, 89, 144,. . . . The ratio between two consecutive Fibonacci 
numbers approaches the golden ratio when these numbers become very large. For small 
lattice sizes it is better to redefine the parameter j3 to be such a ratio of Fibonacci numbers. 
We do this in order to get better results for smaller sizes like L = 8 or 13, which become 
relevant for numerical solutions in the presence of interactions. It is worth mentioning that 



departure from the critical A c = 2 was observed in the experiments of Roati et al. js]. They 
find localization for A > 7, using an incommensurate potential with /3 ~ 1.1972. This value 
of /3, being significantly different from a ratio of consecutive Fibonacci numbers (the golden 
ratio is ~ 1.6180), does not ensure the duality condition [27j that yields A c = 2. 
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FIG. 3: (Color online) Top: Average entanglement as a function of the disorder strength A for 
different lattice sizes in the Anderson model. Averages were calculated from 5000 random con- 
figurations of the potential. Bottom: Derivative of the entanglement with respect to the disorder 
strength A for the same sizes. The location of its minimum for a given system size L is taken as the 
corresponding critical disorder strength A^. On both panels, curves are associated to increasing 
lattice sizes from top to bottom. 

B. Entanglement 

The concept of quantum entanglement can be used to study localization if we choose a 
basis for the Hilbert space composed of states that are localized at each of the lattice sites. 
Then, a fully extended ground-state will have maximum entanglement of the basis states. 
In the opposite limit, the entanglement will be zero for localization at a single site. So, the 
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FIG. 4: (Color online) Entanglement (top) as a function of the potential amplitude A for chosen 
lattice sizes in the AA model, and the corresponding derivative (bottom). Notice the sudden drop 
in the first case, and sharp minima in the second near the point A = 2, except for sizes that are 
not Fibonacci numbers, here exemplified by L = 54. The size effect is better noticed in the bottom 
panel, where deeper minima correspond to larger lattice sizes. 

ground-state Shannon entropy is a good measure of entanglement. It is given by 

S = -^Pilog 2 Pi, (5) 

i 

where pi = |q| 2 , with Cj representin g th e coefficient of the i-th basis state in the expansion 



of the ground-state vector \ty ) [30|, |3l|. For a lattice of L sites, the maximum value of S 
is log 2 L, and occurs when all the Pi are equal. We define S = SJ log 2 L as our measure of 
entanglement, whose maximum value will be S miLX = 1- 

In Fig. [3] (top panel) we show the ground-state entanglement S (averaged over disorder) 
for the Anderson model as a function of the disorder strength A for various lattice sizes 
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L. We define the characteristic Al in this the inflection point of the entanglement 

curve. It can be better visualized as the minimum in the derivative of S with respect to A, 
as shown in the bottom panel of Fig. [3J Notice that the minimum moves toward A = as 
L increases. We will come back to this point in our finite-size-scaling analysis of Sec. IIHI 

In the case of the AA model, the ground-state entanglement, when calculated for lattice 
sizes belonging to the Fibonacci series, shows a sudden drop near A = 2, as can be seen in 
Fig. H] (top panel). This indicates that the derivative dS/dA is again a good marker of the 
quantum phase transition. The bottom panel of Fig. H] shows this derivative as a function 
of A for different lattice sizes. We can see a sharp minimum in the derivative essentially at 
A = 2, more or less independently of the lattice size. Here, the size effect appears mainly 
in the depth of the minimum, which increases dramatically with increasing size. A non 
Fibonacci size has been included for comparison in Fig. HJ 



C. Fidelity 

The fidelity is a measure of how distinguishable two quantum states are. For pure states, 
it is defined as the absolute value of the overlap between these two states. We use here the 
scalar product of two ground-state vectors of Hamiltonians with slightly different values of 
A, writing the fidelity as 

J-(A) = |(tf (A - 5A)|^ (A + 8A))\ , (6) 

which has values between and 1. Generically, the fidelity exhibits a minimum at a critical 
parameter value characterizing a QPT. We have checked that although the choice of SA 
affects the magnitude of the minimum, which can be made arbitrarily small, the value of A 
at which the minimum occurs is independent of 8A. Here we use 8 A = in all numerical 
calculations. In the case of the Anderson model, we perform the configuration average on 
J 7 (A), which means that both ground states are obtained for the same realization of the 
random local energies. 

The average ground-state fidelity for the Anderson model is shown in Fig. |5] (top panel) 
as a function of the disorder strength for different lattice sizes. It is interesting to notice 
that it has a minimum for all lattice sizes at A = 0, which is the expected critical disorder 
strength in one dimension for L — >■ oo. The inset of Fig. [5] shows that this remains true even 
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FIG. 5: (Color online) Average ground-state fidelity (top) in the Anderson model as a function of 
disorder strength for different lattice sizes, and the corresponding derivative (bottom). Averages 
were calculated from 5000 random configurations of the potential. The top-panel inset shows (in 
a much smaller vertical scale) that the minimum at A = exists even for small lattice sizes. 
The position of the maximum derivative for a given system size L is taken as the critical disorder 
strength A^. Deeper minima of the fidelity (higher maxima of its derivative) correspond to larger 
lattice sizes. 

for sizes as small as L = 8 and 10. The finite-size effect is noticeable in the minimum depth, 
which increases with L. Nevertheless, the minima are relatively broad, and an inflection 
point exists, which moves down in A as L increases. This is better seen in the derivative of 
T{ A), shown in the bottom panel of Fig. [5j Similarly to what we observed for the superfluid 
fraction, here the maximum of dJ-'/dA locates a characteristic disorder strength A^ for each 
lattice size. A detailed analysis of the behavior of A^ for large L will be presented in Sec. IIHI 
The fact that ^(A) can signal the critical disorder strength for all lattice sizes is confirmed 
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FIG. 6: (Color online) Fidelity as a function of the disorder strength A for different lattice sizes. 
There is a sharp peak close to the point A = 2. The inset shows two small lattice sizes in a much 
narrower scale close to T = 1. We redefined /3 as 13/8 and 21/13 for L = 8 and 13, respectively. The 
curves indicated at the bottom right show minima that become more pronounced as L increases. 

for the AA model. Figure [6] shows the fidelity as a function of A for different lattice sizes 
(all belonging to the Fibonacci series, except L = 54). A sharp minimum appears very close 
to A = 2 for all the appropriate sizes, clearly indicating the localization transition. In this 
case, one does not extract more information from the derivative. The size effect is evident 
in the minimum depth, which is strongly dependent on the value of L. 

We can observe that the ground-state fidelity is a very precise tool to indicate the local- 
ization transition, even for small lattice sizes. This is specially relevant in the interacting 
case, where numerical solutions must be restricted to relatively small sizes. It is interest- 
ing to mention that the concept of fidelity was also efficiently applied to identify quantum 



phase transitions in other models in the BEC scenario [32|, |33[. For the Bose- Hubbard model 
without disorder, it was shown [20| that the fidelity is a clearer indicator of a QPT than the 
entanglement. 



III. FINITE-SIZE SCALING 



In the previous section, we have seen that good signatures of the localization transition 
in the one-dimensional AA model are provided by the superfluid fraction, the ground-state 
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FIG. 7: (Color online) Logarithmic plot of the critical disorder strength Ax in the Anderson model 
as a function of the lattice size, evaluated from superfluid fraction (squares), derivative of the 
entanglement (circles), and derivative of the fidelity(triangles). The lines are fittings to the scaling 
law (J7J)) which is satisfied with an exponent 7 ~ 1.5 and with A^ = for all three quantities. 

entanglement, and the ground-state fidelity. All these quantities (or appropriate derivatives) 
correctly locate the critical amplitude of the incommensurate potential (equivalent to a 
"disorder strength"), more or less independently of lattice size, provided this size is not too 
small, and is chosen to be nearly commensurate with an integer number of oscillations of 
the AA potential. 

For the true Anderson model, however, the results are less clear-cut upon a simple di- 
rect visualization. We denned characteristic disorder strengths Al in Sec. HIl with different 
definitions depending on the physical quantity evaluated. In order to determine more quan- 
titatively their size dependence, these different A^'s where calculated for lattice sites from 
6 to 300, with random averages being performed typically over 5000 configurations. Despite 
the different definitions, if we make a log- log plot of Al vs L , as is done in Fig. U\ we 
see that in all the three cases (superfluid fraction, entanglement, and fidelity) the following 
scaling law is obeyed: 

A L -A 00 = CL~\ (7) 

where C and 7 are positive real numbers, and Aoo is the critical disorder strength in the 
thermodynamic limit. Despite the fact that the prefactor C in Eq. (J7|) differs for the three 
definitions of A^, we find the same exponent 7 ~ 1.5 in all three cases, and the results are 
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consistent with Aoo = 0, in accordance with the exact value for Anderson localization in 
one dimension. It is also worth noticing in Fig. [7] that the scaling law (J7]) is satisfied even 
for lattice sizes as small as L = 6. 

IV. CONCLUSIONS 

We presented a detailed numerical study of disorder-induced localization in the one- 
dimensional Anderson and Aubry-Andre models of Bose-Einstein condensates. Through 
exact numerical diagonalization of the Hamiltonian matrix, we studied the superfluid frac- 
tion, entanglement, and fidelity as quantities that can signal the localization transition even 
at finite size. Here we restricted our analysis to the non-interacting limit in order to reach 
large lattice sizes. With this, we were able to verify that finite-size-scaling laws are obeyed 
by all the studied quantities in the case of the Anderson model, with a common exponent. 
In contrast, finite-size effects in the Aubry-Andre model are less pronounced, and mainly 
related to the incommensurability between the potential and the system size. If the incom- 
mensurate potential is tuned so that an integer number of periods nearly coincides with the 
system size, then the localization transition is fairly sharp at the critical potential strength. 

Even though all the three quantities analyzed here show clear signatures of the localiza- 
tion transition, we find that the fidelity is less prone to finite-size effects. This makes it 
specially suited to be used in the interacting case, where numerical solutions are restricted, 
for practical reasons, to relatively small sizes. 
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